a rope slides frictionlessly down the edge of a table

[Analog_Engine.git] / scripts / TP1 03.08 Charge in em field.AESL

diff --git a/scripts/TP1 03.08 Charge in em field.AESL b/scripts/TP1 03.08 Charge in em field.AESL
new file mode 100644 (file)
index 0000000..e56c9bd
--- /dev/null
+++ b/scripts/TP1 03.08 Charge in em field.AESL
@@ -0,0 +1,33 @@
+IDENTIFICATION DIVISION
+PROGRAM-ID ChargeInEMfield
+VERSION 20240116
+COMMENT A mass m with charge q>0 is located in a static electro-magnetic field [0,0,E], [0,-B,0].
+COMMENT Find the trajetory
+COMMENT Differential equation (Lorentz force): x'' = omega*z', z'' = a*omega^2 - omega*x'
+COMMENT good start settings are omega=0,92, a*omega^2=0,05
+
+ENVIRONMENT DIVISION
+ENGINE Anabrid-THAT
+TIMEBASE 1ms # corresponds to 1s
+REQUIRES INTEGRATOR 4, INVERTER 1, COEFFICIENT 3
+
+DATA DIVISION
+OUTPUT OUTPUT.X x
+OUTPUT OUTPUT.Z z
+COEFFICIENT.1 OMEGA_X # -omega for x'
+COEFFICIENT.2 A*OMEGA^2 
+COEFFICIENT.3 OMEGA_z # omega for z'
+
+PROGRAM DIVISION
+x'' -> INTEGRATOR -> -x'
+-x' -> INTEGRATOR -> x
+-x' -> COEFFICIENT.OMEGA_X -> -omega*x'
++1  -> COEFFICIENT.A*OMEGA^2 -> a*omega^2
+a*omega^2, -omega*x' -> INTEGRATOR -> -z' # input is z''
+-z' -> INTEGRATOR -> z
+-z' -> INVERTER -> z'
+z'  -> COEFFICIENT.OMEGA_Z -> omega*z' = x''
+
+OPERATION DIVISION
+MODE REPEAT
+OP-TIME 21ms